Coherence measures for heralded single-photon sources 
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Single-photon sources (SPSs) are mainly characterized by the minimum value of their second- 
order coherence function, viz. their g^ ' function. A precise measurement of g^ ' may, however, 
require high time-resolution devices, in whose absence, only time-averaged measurements are acces- 
sible. These time-averaged measures, standing alone, do not carry sufficient information for proper 
characterization of SPSs. Here, we develop a theory, corroborated by an experiment, that allows 
us to scrutinize the coherence properties of heralded SPSs that rely on continuous- wave parametric 
down-conversion. Our proposed measures and analysis enable proper standardization of such SPSs. 
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The demand for ultra-secure communication, high- 

Secision measurement, and super-efficient computation 
has resulted in the emergence of optical sources that 
create single — and, ideally, only single — photons in a 
heralding and/or on demand way [3] . The rapid progress 
in this area has been followed by its early introduction 
to the commercial market |3|, even before finalizing a 
proper set of standards for characterizing such devices. 
For a heralded single-photon source (HSPS) that relies 
on the spontaneous parametric down-conversion (SPDC), 
where the detection of idler photons heralds for the pres- 
ence of signal photons, two figures of merit are generally 
of crucial importance. The first is the ternporal corre- 
lation between the signal and idler beams [J|, on which 
our heralding mechanism relies, and the second is the 
second-order degree of coherence for the heralded signal 
photons [5[ . The challenge of measuring either of these 
figures lies in the large bandwidth of the SPDC process, 
thereby the very narrow widths of such correlation func- 
tions. In fact, what we can commonly measure in an 
experimental setup is only a time-averaged version of the 
actual figure. It is important then to recognize all major 
parameters that affect our measurement results, and put 
them together to come up with well-defined, and read- 
ily measurable, figures of merit for HSPSs. This Report 
carves into the theoretical aspects of such problems and 
addresses the above coherence measures, and their corre- 
sponding time-averaged figures, with an accuracy never 
presented before. Our analysis accounts for the contri- 
bution of multi-photon states in the SPDC process as a 
function of pump power, or, effectively, the single-photon 
generation rate, for different widths of the coincidence 
window and photodetectors' time resolutions. Our the- 
oretical predictions are corroborated by our experimen- 
tal results, and they together provide a prescription for 
proper characterization of HSPSs. 

Figure [T] shows the setup of our HSPS along with the 



Hanbury Brown and Twiss (HBT) interferometer used 
for the measurement of the coherence function. In our ex- 
periment, a 405 nm continuous- wave laser pumped a peri- 
odically poled KTiOP04 (PPKTP) crystal. The crystal 
from Raicol was cut to 10 x 2 x 1 mm'^ for propagation 
along the x axis and poled with a 10 pm period to support 
typc-II SPDC, where signal and idler photons have or- 
thogonal polarizations. By using an oven, frequency de- 
generacy was reached at 39 °C with a stability of ±0.1 °C. 
A polarizing beam splitter (PBS) split the two beams into 
two different spatial modes. A photodetection event on 
the idler beam heralded the presence of one or more pho- 
tons on the signal beam, which went through an HBT in- 
terferometer consisting of a 50:50 beam splitter followed 
by two bandpass filters (810 nm, 10 nm bandwidth) and 
multi-mode fibers that coupled the light to photodetec- 
tors. We employed single-photon counting modules from 
Pcrkin-Elmer with nominal dead-times of 45 ns. Typi- 
cal measured photon count rate for our setup was up to 
800,000 counts/s in the idler channel, with a signal-idler 
coincidence rate amounting up to 10% of that value. 

The detection times for the signal and idler beams were 
recorded by a time-tagging unit from Dotfast Consulting 
with a nominal temporal resolution of 156.25 ps. The 
time-tagging unit streamed the time tags to a computer, 
by which we could record single, double, and triple detec- 
tion events between the three channels, i, s\, and S2, in 
Fig. [l] Coincidence windows were implemented only in 
software. The complete system of photodetectors, power 
supplies, time-stamping electronics, and the USB inter- 
face fit in a 30 X 30 X 30 cm^ box. 

The first coherence measure studied here is the tem- 
poral correlation between signal and idler beams, as a 
measure of reliability of our heralding mechanism, and is 
defined as follows 
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FIG. 1: (Color online) Experimental setup for our heralded 
single-photon source. A blue laser is cleaned from infrared 
fluorescence by a color filter, is focused by cylindrical lenses 
onto a PPKTP crystal (cut and poled for type-II), and is 
then removed from the down-converted beam by a dispersive 
prism. A polarizing beam splitter (PBS) splits the down- 
converted beam into the signal and idler arms. The idler 
beam i is used as a trigger and, in an HBT interferometer, 
the signal beam is split by a 50:50 beam splitter into si and 
s2 for the coherence function measurement. Bandpass filters 
block background light from entering the multi-mode fibers 
that couple the light to the single-photon detectors. 



where Es{t) and Ei{t) represent the scalar photon- units 
positive-frequency field operators for the outgoing signal 
(s) and idler (j) beams, respectively. The joint state of 
signal and idler is a zero-mean Gaussian state whose only 
nonzero second-order moments are given by [6| 

e*"--/2i?(T) = {El{t + T)Ekit)), k = .s,i, (2) 
e-^^v(t+r)/2(j{^r) = {E,{t + T)E,{t)), (3) 

which represent the auto-correlation function and the 
phase-sensitive cross-correlation fimction between signal 
and idler fields, respectively, and where u)p is the pump 
frequency. All other moments can be obtained by us- 
ing the quantum form of the Gaussian moment- factoring 
theorem [7|. For instance, the numerator in Eq. ([T]) can 
be simplified as follows 
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which represents the coincidence rate for having a signal 
photon at time t + t and an idler photon at time t. 

In the low-gain regime of parametric down-conversion, 
which prevails in our case, the auto- and cross-correlation 
functions can be approximated by the following cxprcs- 



R{t) 



i?SPDc(l+T/At) -At<r<0 
i?SPDc(l-T/At) 0<r<At 
elsewhere 
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where i?sPDC is the rate of photon generation for the sig- 
nal/idler beam and 1/At is the bandwidth of the SPDC 
process, and 
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with the difference in the speed of light for ordinary and 
extraordinary axes in the crystal being compensated. 

From Eqs. ^ and (g]), we have gfi{t,T) = 1 + 
\C{t)/R{Q)\^ . For our source, at 50 mW pump power, 
i?(0) = i?sPDC « 43 MHz and 1/At « 3THz, which re- 
sults in a peak value of 1 -f l/(Ati?sPDc) ~ 7 ■ 10* for 

f2l 

9si V',t) at t = 0. The coherence function quickly drops 
to its minimum value one within a sub-picosecond period, 
however the finite time resolution in our experiment will 
smooth this feature out as we show next. 

In order to measure g),^ (t,T), we approximate Psi{T) 
by the rate of coincident events, Nsi{T), in which an idler 
photocount is observed at time t and a signal photocount 
is observed in the interval [t+r — Tcoin, ^+T-f Tcoin], where 
2Tcoin is the width of our chosen coincidence window. 
Because of the photodetectors' time jitters, and neglect- 
ing dark counts throughout the Report, a photodctcction 
event at time t only implies the existence of one or more 
photons in a neighborhood around time t. For simplic- 
ity, we assume that the detection time corresponding to 
a photon that hits the detector's surface at time t is uni- 
formly distributed over the interval [t—Td, t+Td\, where r^ 
is the time resolution of the photodetectors. We can then 
write the observed value for Ngi (r) in terms of Pgi (r) in 
the following way [8| 
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is the coincidence rate for detecting a signal photon at 
time t + T and an idler photon at time i, where u{t) = 
\/{2Td) if \t\ < Td, and zero otherwise. 

Figure[2]shows the experimental and the theoretical re- 
suits for the time-averaged coherence function gl^ (r) = 
Nsi{T)/R'^{0) for different values of Tcoin- Experimen- 
tally, R^ (0) was determined by the product of the signal 
and idler count rates. For the theoretical graphs, we 
used the low-gain correlation functions given by Eqs. ^ 
and (0 with i?sPDC = 43 MHz and 1/At = 3THz. From 
Eqs. (O-®, we see that Nsi{T) has an almost fixed value 
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FIG. 2: (Color online) Top: Measurements (symbols) and 
theory predictions (lines) of the time-averaged coherence func- 
tion gli (t) for the signal and idler photons. The low-gain 
regime theory curves are in striking agreement with the data 
using the following parameter values i?sPDC ~ 43 MHz, r^ — 
350 ps and 1/At = 3THz, where the last one was measured 
separately by spectroscopy and the other two were adapted 
to provide the subjective best visual fit simultaneously to all 
three sets of data. Bottom: Same data as top but with 
magnified ordinate in order to reveal structure that is caused 
by double refiection of photons from the fiber input endface 
and the surface of another optical element in the experiment. 
While this effect causes only very small deviations (0.2% of 
the central peak) from the expected flat line, these deviations 
cause noticeable ringing of the measured gc (see Fig. |3]). 



for T e [-Tcoin + At -I- Td, Tcoin - At - Trf], inversely pro- 
portional to Tcoin- As we get farther from the center, the 
time-averaged coherence function drops to its minimum 
value one as expected. The theoretical graphs are in 
striking agreement with our experimental results, which 
clearly demonstrate the strong temporal correlation be- 
tween signal and idler beams. 

The second coherence measure that we consider here is 
the degree of second-order coherence for the signal field, 
conditioned on observing an idler photocount at time ti, 
defined as follows 



gP{h,h\U) = 
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where (•)pm is the average over the post-measurement 
state assuming sufficiently high time resolution and unity 
quantum efficiency for the idler photodetector. 

To model the measurement on the idler field opera- 
tor, we use a heuristic approach in which a photode- 
tection event at time U on the idler beam is modeled 
by the measurement operator, [9|], Ei(ti). We can show 
that if we allow for infinitely high time resolutions this 
method provides us with the correct result [lO|. The 
post-measurement averaging, for any operator X, will 



then be given by 

(X)p,„ = {e}{U)XE,{U))/{e}{U)E,{U)). (10) 

The conditional coherence function in Eq. ([9]) can then 
be written as follows 

Psi{tl — ti)Psi{t2 — ti) 

where, using the quantum version of the Gaussian 
moment-factoring theorem along with Eqs. ([2]) and ([3]), 

P^f\ti,t2,U) = {EJiU)Elih)Elit2)E,it2)E,{h)E,iU)) 

= i?(0)[i?2(0) + |i?(ti-f2)P] 

+ 2di{C{ti-t,)C*{t2-U)R{ti-t2)} 
+ i?(0)[|C(ii-tOP + |C(t2-f,)l'] (12) 

is the multi-coincidence rate for having signal photons at 
times ti and i2 and an idler photon at ti. 

The first of several interesting special cases we consider 
is the coherence function at the trigger time, i.e., 
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It is clear that if R^{0) <C \C{0)\^ then gr'iti,t,\t,) « 0. 
This is the same requirement that we had for observing a 
large g^j (t, 0) in Eq. ([T|). and, therefore, a low value for 
9c {ti,ti\ti) is guaranteed if 3^^ (i,0) » 1. The second 
interesting case is when ii = ti but \t2 — ti\ ^ 2 At. In 
this case, gc (ti,t2\ti) ~ 1 provided that R^{t2 — ti) « 
and \C'^it2 - tj)| « 0. This implies that our HSPS has a 
coherence time on the order of At. Finally, let us consider 
the case when \ti — ti = t2 — ti\ ^ 2Af, i.e, when there 
is no correlation between the trigger time and the signal 
beam. In this case, gc {ti,t2\ti) « 2, which is expected 
because, in the lack of any triggering event, both sig- 
nal and idler beams individually obey the thermal-state 
statistics, for which the second-order coherence function 
has a maximum value of two [lll | . 

To quantitatively characterize our HSPS, it is interest- 
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For an ideal HSPS, we expect that gc (0) = 0. In our 
case, from Eq. HH), gi^\o) « 6 ■ 10"^ < 1. Here, we 
measure a time-averaged version of the coherence func- 
tion by approximating Psiir) with Ngilr) as before and 
Pgj (0, r, 0) with iVg^j (r), the count rate for triple coin- 
cidences of an idler photodetection event at t^ = 0, and 
two signal photodetection events at ii G [—Tcoin, '''coin] 
and t2 € [t — Tcoin, t + Tcoin]- In our HBT interferom- 
eter, we can equivalently look at the number of triple 
coincidences on the idler and S1-S2 photodetectors. By 
accounting for the resolution of the three photodetectors 
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FIG. 3: (Color online) Measured (symbols) and calculated 
(lines) time-averaged conditional coherence function cjc (t). 
The theory lines were calculated using the same parameter 
values as in Fig. [2] The purely statistical errors of our data 
are on the order of the symbol size in the figure and therefore 
not shown. As explained in the caption of Fig.[2]photons that 
are refiected twice cause the apparent ringing. 



involved in our measurement, we obtain 
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is the multi-coincidence rate for detecting an idler photon 
at time and two signal photons at times ti and t2- 
Figure [3] shows our measurement results for the 

(2) 

time-averaged conditional coherence function cjc (t) = 
N^f\T)R{0)/[N,,{0)N,,{T)] for three different coinci- 
dence windows, which result in three different widths for 
the observed central dip. Here, R{0) is obtained from 
the idler count rate in the experiment. As explained in 
Fig. [21 the ringing structure in Fig.[3]is caused by double 
optical reflections. The graphs, nevertheless, exhibit the 
signature of a good SPS as the measured value of ^c (0), 
at 43 MHz single-photon generation rate, in Fig. [3l is 
0.0781 ± 0.0006 for 2rcoin = 0.78 ns and Td = 0.35 ns. 
By reducing the pump power we can reduce this resid- 
ual 5c (0) almost arbitrarily at the expense of reducing 
the total count rate. To see how the depth of the dip in 
Fig. Ovaries with the coincidence window, in Fig. [4l we 
have plotted gc (0) versus 2Tcoin- It can be seen that, 
for Tcoin ^ Td, (jc (0) is determined by t^, whereas, for 
Tcoin ^ Td, it is almost linearly increasing with Tcoin- 
Our theoretical treatment is again well capable of repro- 
ducing the measurement results. The graph shown in 
Fig. m exemplifies the fact that a single value for (jc (0) 
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FIG. 4: (Color online) Experimental (symbols) and theoreti- 
cal (line) results for the minimum of the time-averaged con- 
ditional coherence function, g^ (0), as a function of the coin- 
cidence window 2rcoin using the same set of parameters as in 
Fig. [5] The dashed line is for ideal photodetectors (r^ = 0). 



does not bear enough information to quantify the source 
performance. At a fixed rate, the interplay between the 
coincidence window and the time resolution of photode- 
tectors must also be accounted to give a proper figure of 
merit for an SPS. 

In this Report, we theoretically and experimentally 
studied the coherence properties of heralded single- 
photon sources that use parametric down-conversion. 
We used the Gaussian characteristics of down-converted 
fields to analytically find the degree of second-order co- 
herence between signal and idler fields as well as for the 
signal field, individually, when it is conditioned on the 
detection of an idler photon. Our theory is well capable 
of reproducing our experimental results, which demon- 
strated a high-quality source of sub-picosecond single 
photons. It also allowed us to study the impacts of the 
chosen coincidence window, the down-conversion param- 
eters, and the resolution of photodetectors on the out- 
come. Such analysis enables proper standardization of 
such devices. We would like to thank N. Liitkenhaus, I. 
SoUner, and A. Safavi-Naeni for their technical assistance 
and acknowledge NSERC, CFI, ORF-RI, ORDCF, Quan- 
tum Works, CIPI; and CIFAR for their financial support. 
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